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The nonlinear magnetic X3 an< i spin-glass \sg susceptibilities in zero applied field are obtained, 
from tempered Monte Carlo simulations, for three different spin glasses (SGs) of Ising spins with 
quenched site disorder. We find that the relation — T 3 X3 = Xsg — 2/3 (T is the temperature), 
which holds for Edwards- Anderson SGs, is approximately fulfilled in canonical-like SGs. For nearest 
neighbor antiferromagnetic interactions, on a 0.4 fraction of all sites in fee lattices, as well as for 
spatially disordered Ising dipolar (DID) systems, — T S X3 and Xsg appear to diverge in the same 
manner at the critical temperature T Bg . However, — T 3 \3 is smaller than Xsg by over two orders of 
magnitude in the diluted fee system. In DID systems, —T 3 X3/Xsg is very sensitive to the systems 
aspect ratio. Whereas near T sg , \ sg varies by approximately a factor of 2 as system shape varies 



from cubic to long-thin-needle shapes, \3 sweeps over some four decades. 
PACS numbers: 75.10.Nr, 75.50.Lk, 75.30.Kz, 75.40.Mg 



I. INTRODUCTION 



The existence of an equilibrium phase transition into 
the spin glass (SG) phase has not yet been convincingly 
established for some spin glasses. The development of 
the parallel tempered Monte Carlo (TMC) algorithm 1 
has enabled one to observe, bypassing anomalously long 
relaxation processes, SG models in equilibrium at low 
temperatures. Thus, correlation lengths £ have been de- 
termined from the equilibrium behavior of {{siSj)^) q , 
where si = ±1 is for a spin at site i, and (• ■ • )t and 
(• • • ) q stand for a thermal average and for an average 
over quenched randomness, respectively. There is evi- 
dence, from Monte Carlo simulations, that £ grows as 
linear system size L in (i) the Edwards-Anderson (EA) 
model 2 ^ at some nonzero temperature T sg in three di- 
mensions, in (ii) geometrically frustrated systems, such 
as strongly site-diluted Ising models, with nearest neigh- 
bor antiferromagnetic (AF) bonds, on foe lattices^ and 
in (iii) strongly site-diluted Ising models with dipole- 
dipole interactions, such as in LiHo^Yi-zF^ We refer 
to the latter systems as disordered Ising dipolar (DID) 
systems!^ At least for DID systems, some numerical ev- 
idence that is unfavorable for the existence of a phase 
transition also exists^ The divergence of £ implies the 
divergence of the so called spin-glass susceptibility Xsg 
at T sg , where \s g = N^ 1 J^ij ((siSj^g an d is the 
number of spins. 

Convincing experimental evidence for the existence of 
an equilibrium phase transition into the SG phase is 
harder to obtain. This is mainly (i) because very long 
relaxation processes make equilibrium observations diffi- 
cult, and (ii) because neither £ nor \sg can be directly 
observed. Instead, the SG transition is usually charac- 
terized by the nonlinear magnetic susceptibility X3^ It 
is defined by 



assuming m(—H) = —m{H). Canella and Mydoshii 
were first able to measure (in gold-iron alloys) huge values 
of X3'- Ts g x?,/xi ~ 10 5 (from here on, we let Boltzmann's 
constant and Bohr's magneton equal 1) near T = T sg . 
Later, X3 was showni^ to diverge in other canonical SGs 
as a power of T — T sg . For the EA model^ origi- 
nally inspired by the discovery of Canella and Mydosh, 
ChalupaM showed long ago that 



X3 = T-% X sg - 2/3) 



(2) 



if no field is applied. Thus, the critical behavior of X3 and 
Xsg, which one observes in simulations, can, at least for 
the EA model, be clearly related to the critical behavior 
of X3j which one observes experimentally. 

The three models we study are governed by the Hamil- 
tonian, 



U 



ijX%SiX j Sj , 



(3) 



XiH + X3 H 3 



(1) 



where the sum is over all i and j lattice sites, 3\j is model 
specific, Xi = ±1 is a quenched random variable, and Si 
is a (±1) Ising spin at site i. All sites are occupied with 
the same probability, x = (xi) q , where the q subscript 
stands for a quenched average over all site occupancy 
arrangements. 

The aim of this paper is to find how — X3, an experi- 
mentally measured quantity, and Xsg, a quantity which 
is often calculated, are related in site-diluted SGs. More 
specifically, numerical results from TMC are sought for 
(i) Ising spins, with Ruderman-Kittcl-Kasuya-Yoshida 
(RKKY) interactions^ which are randomly located on 
a small fraction of all lattice sites, (ii) a geometrically- 
frustrated Ising spin system, mainly, randomly located 
Ising spins, with nearest neighbor AF interactions, on a 
0.4 fraction of all sites of an fee lattice, and (iii) DID 
systems on a small fraction of all lattice sites. The out- 
come of these calculations is unknown, because Eq. @ 



2 



has not been derived for site diluted SGs. On the other 
hand, X3 an d Xsg can exhibit the same critical behavior 
in site diluted SGs if they and the EA model belong to 
the same universality class. This has been predicted 1 ^ to 
be so for the first of the above three models, but not so, 
as far as we know, for the other two models i 17 ' 18 

An outline of the paper follows. Details about the 
procedure we follow in our calculations are given in Sec. 
IITI For various sizes of each of these systems, we ob- 
tain £/L, Xsg, and \3- Data for £/L is used to establish 
the phase transition temperature T sg between the para- 
magnetic and SG phases. We then compare how Xsg 
and X3 vary with system size and with temperature in 
the vicinity of T sg . The results obtained for each sys- 
tem are given in each of the subsections of Sec. IIIII 
Very briefly, these results follow. We find that \3 a P~ 
proximately follows Eq. (TTJ) in strongly diluted systems 
of Ising spins with RKKY interactions. On the other 
hand, — T 3 X3 is a over a couple of orders of magnitude 
smaller than Xsg in a (x = 0.4) site-diluted AF Ising 
model on an fee lattice. Nevertheless, both X3 an d Xsg 
appear to have the same critical behavior. Finally, in 
DID systems, — T 3 X3 an< A Xsg seem to diverge similarly at 
T sg . However, —T 3 X3/Xsg varies sharply with systems' 
shape. Taking into account demagnetization effects, we 
estimate in Sec. HVI how — T 3 X3 varies with system shape 
for high aspect ratios. Near the transition temperature, 
— T 3 x^/Xsg increases from — T Xz/Xsg ~ 10 -2 for cubic 
shape systems to — T 3 X3/Xsg ~ 10 2 for very thin long 
prisms. Our conclusions are summarized in Sec. [V] 

As a byproduct, we have obtained values for 77 and T sg 
in these three systems. They are listed in Table II. 

From here on, in addition to ks = 1, (J-b = 1, we let 
m = N~ z2i((si)T)q) an d assume spins in all models 
point up or down along the z-axis, sometimes referred to 
as the magnetization axis. 



II. PROCEDURE 



pies, and noting (i) that both (M 4 ) and (M 2 ) 2 involve 
sums over four-spin terms, such as J2ijki(( SiS j SkSl ) T )q 
an d J2t 3 ki{{ s ^ s j)T{skSi}T)q, respectively, and (ii) that 
any term in which one or more subindices is unpaired 
vanishes. To see this, assume the k index is unpaired in 
either of the two sums over ijkl indices. Now, consider 
all exchange bonds Jkm between the fc-th and any other 
site. Let's assume the probabilities for J^ m and —Jkm 
for all m while all other exchange constants remain un- 
changed are equal. (This requires exchange bonds to be 
independently random.) It follows that the probabilities 
for Sk and — Sfc, for any given configuration of all the 
other spins, are equal. This is the gist of the proof. For 
further details, see Ref. 14 1. The proof fails for site- 



diluted systems, because exchange bonds are not then 
independently random. 

We simulate a set of identical systems at tempera- 
tures T min , T min + AT, T min + 2 AT, ...T max follow- 
ing standard TMC rulesJ* We choose T max ~ 2.5T sg , 
Tm m ~ 0.5T sg , and all AT such that at least 30% of all 
attempted exchanges between systems at T and T + AT 
arc successful for all T. We let each system equilibrate for 
a time t s and take averages over an equally long subse- 
quent time t s . Time t s satisfies two requirements, (i) that 
((M) 2 T ) q «0.1((A/ 2 ) T } 9 obtains for all Te [T min ,T max ], 
and (ii) that, systems which start from either random 
configurations or from (assumed) equilibrium configura- 
tions come to the same condition, as specified in Ref. 
[l9l |. after time r s . Values of r s , of the number of sam- 
ples N s with different quenched randomness over which 
averages were taken, and of the site occupancy rate x, 
are given in Table I. 

Periodic boundary conditions are used throughout. 
For DID systems, we make use of Ewald sums^S 



For the correlation length £, we make use, as 
is become 
definition^ 2 - 



has become standard practice ; 2 ' 8 ' 19 ' 21 of the original 



To calculate X3, we make use of 

6X3-^ 1 T- 3 ((A./ 4 ) T -3(M 2 ) 2 ,), (4) 

where M = Nm, which holds for H = 0. This equa- 
tion follows from Eq. (TTJ) by (i) repeated differentia- 
tion with respect to H of the canonical ensemble average 
expression for and by (ii) letting A = 0, where 

A = -4(M 3 )(M) + 12(M 2 )(M) 2 - 6(M) 4 . This is justi- 
fied for finite systems with up-down symmetry if averages 
are taken over infinite times, since (M n ) = then for all 
odd n. The order in which system sizes and averaging 
times are taken to infinity is irrelevant for the param- 
agnetic phase. Equations @ and (0]), as well as all the 
results below are only claimed to hold for T > T sg . Note 
that Eq. ([J) is valid for each realization of quenched 
disorder. Chalupa derived Eq. © from Eq. Q for 
the EA model by first averaging over all system sam- 



1 



4sin 2 (fc/2) 



Xsg 



Xs 9 (k) 



f 



(5) 



where Xs 3 ( k ) = N 1 J2ij(( s i s ])T)<i exp(ik • r^), and we 
let k = (2?r/L,0,0). Note that Xs 9 (k = 0) = Xsg and 
that, as L — > 00, £/L vanishes in the paramagnetic phase, 
remains finite at T = T sg , and either grows without 
bounds below T sg , as in a conventional phase transition, 
or remains finite as in the XY model in two dimensions. 
The point where £/L curves for various system sizes meet 
as T decreases defines T sg for us. 

All systems we report on below have a common feature: 
fractional errors for X3 ar e an order of magnitude larger 
than the ones for Xsg a nd for xi- Unless otherwise stated, 
error bars are given for X3 an d related quantities, but not 
for Xsg or xi, which arc all smaller than icons for their 
data points. 
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TABLE I. The number of samples N s and the number of 
Monte Carlo sweeps r a that were taken both for equilibration 
and for the subsequent averaging times are given in thousands 
of Monte Carlo sweeps. There are L 3 lattice sites in RKKY 
and FCC systems, but LxLx2L lattice sites for DID systems. 
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FIG. 1. (Color online) (a) Semilog plots of £/L vs T for 
[m = ±1) Ising spins with RKKY interactions, randomly 
located, on a 0.1 fraction of all L z sites. The numbers in the 
box are L values, (b) Same as in (a) but for Xsg/L 2 ~ ri vs T, 
for r\ = —0.5. 



III. RESULTS 



All results in this section follow from TMC simulations. 



A. Spatially disordered Ising spins with RKKY 
interactions 

The Hamiltonian is given by Eq. ([3]), with = 
e c (coskrij)(a/rij) 3 , ie, an RKKY— interaction, as in a 
canonical spin glass. We let ka — 2ir, where a is a nearest 
neighbor distance, and e c is an energy in terms of which 
all temperatures arc given in this subsection. We let each 
site is be occupied with x = 0.1 probability. 

Plots of £/L vs T for various system sizes are shown in 
Fig. QJi. Not all pairs of curves cross at the same point. 
Let Tij be the temperature where curves for lengths Li 

= 8 and L 3 = 12. Plots 
0.23, and T 2 , 3 = 0.20) 




FIG. 2. (Color online) (a) Semilog plots of -T^xa/i " vs 
T, for r\ = —0.5, for (m = ±1) Ising spins with RKKY inter- 
actions, randomly located, on a 0.1 fraction of L 3 sites, for 
the values of L, which are shown in the box. (b) Same as in 
(a) but for -T 3 x 3 /Xs 9 vs T. 



bounds with system size at T = T sg . Indeed, we note 
that Xsg ~ L 2 ~ v , where rj ~ —0.5 at T sg . This r] value at 
the boundary of the range —0.5 < rj < —0.2 of quoted^ 
values, from Monte Carlo simulations, for the EA model. 
We note in passing that this model and the EA model 
have been predicted^ to be in the same universality class. 

How — T 3 \3 behaves near T sg is shown in Fig. (2)}. It 
varies with T and with L much as Xsg does in Fig. [TJd. 
The plots shown in Fig. [^b are consistent with — T 3 X3 ~ 
Xsg- (The value rj = —0.5 follows from the plots shown 
in Fig. [TJd, not from any fitting of X3 to any desired 
behavior.) 

More significantly, —T^xslXsg appears to approach a 
smooth function of temperature in the neighborhood of 
T — T sg as L — > oo. This is the basis for the main 
conclusion of this section, namely, that — T 3 X3 an d Xsg 
have the same critical behavior. 



B. Site-diluted AF Ising model on a fee lattice 

Each site of a fee lattice is occupied with a (±1) Ising 
spin with a 0.4 probability. The Hamiltonian is given 
by Eq. <j3j> , with = — J if i and j are nearest neigh- 
bors but Jij = otherwise. A 0.4 occupancy rate is 
roughly midway between the lowest value x = 0.195 
for percolation^ in fee lattices and the transition point, 
x ~ 0.75, between SG and AF phases A All temperatures 
are given in terms of J. 

Monte Carlo results for this model arc shown in Figs. 
[3ja, and [3b. We note in Fig. [3ja, that the crossing point be- 
tween pairs of £/L curves drifts leftward as their L- values 



and Lj cross, where L\ = 6, Li 
of T i:j (where T h2 = 0.29, r 1)3 = 

vs 1/LiLj fall onto a straight line which extrapolates to 

T = 0.10 as 1/LiLj — > 0. We thus estimate T sg = 0.10, increase. As for Fig. QJi, let T i} j be the temperature 
and therefore expect T/x = 1.0 for x < 1 values, since 
the 1/r 3 dependence of the interaction implies^ £(x, T) — 
Z(T/x) : 

In Fig. [TJd we can see that Xsg seems to grow without 



where curves for lengths Li and Lj 



and Lj cross, where L\ = 4, 
L 2 = 6, L 3 = 8 and i 4 = 12. A second degree polyno- 
mial fit to a plot of Tij (where T lj2 = 0.70, T lj3 = 0.67, 
Ti, 4 = 0.62, T 2 ,3 = 0.62, T 2A = 0.57, and T 3A = 0.53) vs 
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FIG. 3. (Color online) (a) Plots of £/L vs T for a (x = 0.4) 
site-diluted AF Ising model on an fee lattice of L x L x L sites. 
The numbers in the box are L values, (b) Same as in (a) but 
for Xsg/L 2 -' 1 vs T, for 77 = -0.5. 



0.003 




FIG. 4. (Color online) (a) Plots of -I^xs/^ 2- " vs T, for 
7] — —0.5, for a (x = 0.4) site-diluted AF Ising model on an 
fee lattice of L x L x L sites. The numbers in the box are L 
values, (b) Same as in (a) but for ~T A X 3/Xsg vs T. 



1/LiLj gives Tij — > 0.4 as 1/LiLj — >• 0. We thus estimate 
T sg = 0.4(1), in agreement, within errors, with the values 
found for x = 0.4 in Ref. Q. Plots of Xsg/L 2 '* 1 vs T are 
shown in Fig. (3Jd for 77 = —0.5. This is the best value of 
77 to have Xsg/L 2 ~ v curves for various values of L cross 
at T sg . This value of 77 is, within errors, in agreement 
with the value found for x = 0.4 in Ref. j4|. 

Plots of -T 3 X 3/L 2 ~ n vs T, with 77 = -0.5, are shown 
in Fig. |4ji. The 77 = —0.5 value is taken from the plots of 
Xs g/L 2 ~ ri vs T, not from any fitting of X3 to any desired 
behavior. The curves in Figs. [3)} and [4^, are somewhat 
different, but all curves for L = 6, 8 and 12 in both figures 
within errors, at the same temperature, T sg = 



HJd we notice that — Tx3 <C Xsg, which differs 



do cross, 
0.4. 
In Fig 

markedly from what might have been expected from the 
behavior of the EA model (and from the above results 
for SGs with RKKY interactions). More significantly, we 
observe —Txz/Xsg is, within errors, independent of L for 
the largest values of L, and appears to go into a smooth 



TABLE II. Values for T sg follow from crossing (or merging) 
points of £/L curves for various values of system linear size, 
L. Values of 77 are assigned so that Xsg/L 2 ~ v curves for var- 
ious L values cross at T 3g . Errors in 77 follow from errors in 
T 3g . As explained in Ref. [§], T sg can be obtained for DID 
systems (for all x < 0.5 in sc lattices^ and x < 0.25 ir* 2 ^ 
LiHo^Yi-zF-ft) from the T ag value given below, making use of 
T sg oc x. Similarly for RKKY interactions and all x < 0.1 
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function of T, near T sg , as L — > 00. This suggests that, in 
the thermodynamic limit, both quantities have the same 
critical behavior. 



C. Spatially disordered (±1) Ising dipoles 

Here we consider disordered Ising dipolar (DID) sys- 
tems in sc lattices. We let each site be occupied, with a 
0.35 probability, by a (±1) spin. All spins point up and 
down, along the z-axis. The Hamiltonian is given by Eq. 
©, with 



J v - h dA- I 3-j- - 1 



(G) 



where is the distance between i and j sites, Zij is the z 
component of r.y , hd is an energy, and a is the SC lattice 
constant. 

Let's first recall that, despite some earlier numerical 
evidence to the contrary^ more recent calculations point 
to the existence of a phase transition between the param- 
agnetic and SG phases in diluted Ising dipolar systems^ 
at T sg /x ~ 1 for all x < 0.5. In addition,^ Xsg ~ L 2 ~'' 
at T = T sgi where 77 ~ 0. 

We deal with the magnetic susceptibility here which, 
as is well known, depends on the shape of the system^ 
For this reason we study numerically L x L x nL shaped 
prism systems for various values of n, that is, square-base 
prisms with al:n aspect ratios. 

Plots of £/L vs T/x are shown in Fig. [5^ for n = 2. 
Curves for three different values of L are observed to 
cross at T/x ~ 1.0. This transition temperature value is 
in agreement with the result found in Ref. [U for n = 1, 
mainly, that T sg /x ~ 1.0 for all x < 0.5. 

Plots of Xsg/L 2 vs T/x are shown in Fig. [5Jd. All 
curves are observed to cross at T/x = 0.95. This is ap- 
proximately as in Ref. 0. 

We now turn our attention to the magnetic susceptibil- 
ity. Plots of — T 3 X3/L 2 vs T/x are shown in Fig. [6^ for 
systems of various sizes. These plots resemble the ones 
for Xsg in Fig. [5Jd, but note that the crossing points are 
not quite at the same temperature in the two figures. 

In order to better compare — T 3 X3 an d Xsg, we plot 
in Fig. [5Jd the ratio —T 3 X3/Xsg vs T/x for systems of 
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FIG. 5. (Color online) (a) Semilog plots of f/L vs T/x for 
DID systems on a 0.35 fraction of all L x L x L z sites, where 
1/2 = 2L. The numbers in the box are L values, (b) Same 
as in (a) but for Xsg/L 2 vs T/x. Error bars are smaller than 
icons for all data points. 




FIG. 6. (Color online) (a) Plots of -T 3 xa/L 2 vs T/x for 
DID systems on a 0.35 fraction of all L x L x L z sites, where 
L z = 2L. (b) Plots of -T 3 X 3/Xsg vs T/x for three 1 : n 
aspect ratios, n — 4 (EB and E3) for the two top curves, n = 2 
(full icons) for the three curves in the middle, and n = 1 
(empty icons) for the three lower curves. Error bars hardly 
protrude from icons. For both (a) and (b), the numbers in 
the box are the values of L. 



various sizes with 1 : 4, 1 : 2 and 1 : 1 aspect ratios. For 
a 1 : 4 aspect ratio, only data points for systems with 
4 x 4 x 16 and 6 x 6 x 24 sites appear in Fig. [SJa. A larger 
system with the same aspect ratio would have taken a 
prohibitively long computer time to run. Let L+ be a 
system length such that —T 3 X3/Xsg is approximately size 
independent if L > L+. Clearly, L+ ~ 4 and 6 for n = 4 
and 2, respectively, in Fig. [6}d. For n = 1, ~ 8 seems 
likely. This would be in accordance with the expectation 
that — T 3 X3 an d Xsg have the same critical behavior in 
DID systems, independently of aspect ratio. 

Questions about the sharp variation of —T 3 xz/Xsg 
with respect to aspect ratio naturally arise. What is 
the asymptotic behavior of —Txs/Xsgt This is hard 
to foresee from the data plots shown in Fig. [5Jd. To pro- 



ceed much further numerically is impractical. The next 
section is devoted to this question. 



IV. VARIATION OF X a WITH ASPECT RATIO 
IN DID SYSTEMS 



In this section we derive an approximate equation for 
the variation of X3 with shape in DID systems. 

Consider two systems of the same shape and size. In 
system /, all spin pairs interact. In the other system, 
system t, dipolc-dipole interactions are truncated. In t, 
each spin interacts only with spins that lie within a long 
thin cylinder centered on it, whose axis is parallel to the 
system's z-axis. The radius of the cylinder need not be 
more than a couple of nearest neighbor distances, but 
its length must be much longer than its radius. Let's 
furthermore assume that both systems are homogeneous, 
that is, all sites are occupied (x = 1). Now, we know from 
Ref. [26j that if both systems are in thermal equilibrium, 
and external magnetic fields H t and Hf are applied to 
systems t and /, respectively, such that m is the same in 
both systems, then 



Hf — H t - A„m, 



(7) 



where n comes from / system's 1 : n aspect ratio. Equa- 
tion (JJj) holds because the only effect of the imtruncated 
portion of all dipole-dipole interactions in / is to give 
the so called demagnetizing field, — \ n m. For dipolar 
prisms of 1 : n aspect ratio, a scaling expression, such 
as2£ m = t~P f{Ht~P s ) must therefore be replaced by 

m = r p f[{H-Km)t- ps ]. (8) 

where Hf has been replaced by H . 

Taking the derivative of Eq. ([H]) with respect to H [or, 
more simply, of Eq. ([7]) with respect to m] gives 



1 



Xi(n) Xi(cc) 



+ A, 



(9) 



where %i (n) is the linear susceptibility of a prism with a 
1 : n aspect ratio, and, clearly, dm/dH t = xi(oo). This 
is the well known equation^ that experimentalists^^ 
often use in order to do away with demagnetization ef- 
fects, and thus relate Xi( n )j the measured susceptibility, 
to Xl(oo). 

Taking the d/dH t derivative of Eq. © gives 



1 + A, 



dm 
dHt 



1 dxi(n) 



1 dxi(oo) 



X{(n) dH n xf(oo) dH t 



(10) 



where we have used dH n /dH t = 1 + X n dm/dH t , which 
follows from Eq. (J7J). We next (i) take the d/dH t deriva- 
tive of the above equation, (ii) let H t = H n = 0, and 
(iii) let dxi(oo)/dH t = = dxi{n) / dH n . by up-down 
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-T 3 X /L 2 & y/L 2 
10 2 




FIG. 7. (Color online) (a) Plots of Xi vs l/ n f° r DID systems 
at T = T sg , on a 0.35 fraction of all L x L x nL sites. The 
shown numbers are L values for data points from Monte Carlo 
calculations. The dashed lines follow from Eq. ([9]), assuming 
the three values for xi(oo) that are shown in the box. (b) 
Same as in (a) but for (full icons) —T 3 X3/L 2 and (open icons) 
Xsg/L 2 . The dashed lines follow from Eqs. (J9| and and 
the shown pairs of values, such as 9 and 100, are for Xi(°°) 
and X3(oo) 5 respectively. 



TABLE III. A n values, in terms of hd, for some n in the [0.5, 71 



range. For n > 8, A n 

An ^ xX n - 



For an x site occupancy rate, 



0.5 1 1.5 2 2.5 3 4 5 6 7 
7.419 4.189 2.503 1.611 1.107 0.802 0.471 0.308 0.217 0.160 



symmetry. The result is easily cast into 

X3{oo) 



X3(n) 



[l + A„ X i(ao)] 4 ' 



(11) 



which is the desired expression relating X3( n ) an( i X3(°°)- 
Equations © and (fTTjl enable us to calculate how 
X3i n ) varies with n if we know Xi an d X3 f° r a t least 
an aspect ratio each, as well as A„ . A list of (easily com- 
puted) A„ values for several values of n € [0.5, 7], as well 
as a functional relation for all n > 8, are given in Table 
III. Since the only effect of the long-range portion of all 
dipolc-dipolc interactions in / is to give the demagnetiz- 
ing field^ — A„m, we can calculate all A„ assuming a 
fully occupied lattice in which all spins point up. Since 
only a single state comes into the calculation, no Monte 
Carlo simulation is necessary. This enables us to calcu- 
late X n for very large systems. The fact that only an x 
fraction of lattice sites are occupied in site diluted SGs 
is approximately taken into account by letting A„ — > xX n 
everywhere. Inhomogencities in SGs are thus neglected. 

Equation © gives the three dashed lines shown in Fig. 
[7_k for Xi(°°) = 7, 9 and 11. With two of these values, we 
obtain from Eq. (fTTjl the three curves for X3( n ) shown 
in Fig. [7b for values of X3(°°)- These curves do not fit 
the data points too well. On the other hand, a good fit 



for small system sizes should not be expected. We can 
nevertheless conclude with some confidence that X3( n ) 
does not diverge as n — > oo. Indeed, X3(°°) is most likely 
within the (50, 120) range. Furthermore, observation of 
Fig. [Jjj indicates that X3 {n) / Xsg at T = T sg varies over 
three or four orders of magnitude as system shape varies 
from cubic to infinitely thin needle-like. 
It is perhaps worth pointing out that 



x \ oo 



(12) 



follows immediately from Eqs. © and (fTTjl after Eq. © 
is cast into Xi{ n ) = Xi( 00 )/[l + ^riXi( 00 )]- Equation (fl~2"l) 
implies that X3 sweeps over four times as many decades 
as xi does (compare Figs. [7k and [7b) as n varies. 

Finally, note that the classical or quantum nature of 
DID systems does not play any role in this section. It 
does not matter either whether a transverse field is ap- 
plied, because it does not affected up-down symmetry. 
These equations can therefore be applied, as an illustra- 
tion, to Lii — ccHo^Y4, under a transverse field, as in Ref. 
[24j . where T 2 xs/xi ~ 1 was observed on a 1.6 x 16 x 5 
mm 3 sample. Values of xi a nd X3 that would be some 3 
and 100 times larger, respectively, for a long thin needle- 
like sample can be read off from Figs. [7k. and [7b. 



V. CONCLUSIONS 

By the tempered Monte Carlo method^ we have tested 
whether the relation — T 3 X3 = Xsg — 2/3, which is 
knowni^ to hold for the Edwards- Anderson model, also 
holds for several site-diluted spin glasses of (±1) Ising 
spins, with (i) RKKY interactions, (ii) antiferromagnetic 
interactions between nearest neighbor spins on fee lat- 
tices, and (iii) dipole-dipole interactions. As a byprod- 
uct, we have obtained the values of r\ and T sg that are 
listed in Table II. 

We have found — T 3 X3 ~ Xsg to hold, for Ising spins 
with RKKY interactions occupying a 0.1 fraction of all 
lattice sites. More significantly, —T 3 X3/Xsg appears to 
be (i) independent of linear system size, within errors, 
and (ii) a smooth function of temperature near T sg . This 
suggests — T 3 X3 and Xsg have the same critical behavior. 
Since the RKKY interaction decays as the inverse of the 
cube of the distance, these results must hold for lower 
values of x if the temperature is scaled with x. 

We have found — T X3 to be over two orders of magni- 
tude smaller than Xsg for Ising spins, with antiferromag- 
netic interactions, on a (x = 0.4) site diluted fee lattice. 
Our results are, however, consistent with identical critical 
behavior of these two quantities. 

In DID systems the TMC data (see Fig. [BJd) are con- 
sistent with X3 a nd Xsg diverging in the same manner as 
T — > T sg from above. The sharp variation of —T xs/Xsg 
with aspect ratio, which can be observed in Fig. [SJs, is 
noteworthy. How this comes about from demagnetization 
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effects is explained in Sec. IIV1 In it, relations are derived 
which together with data points coming from TMC sim- 
ulations (see Fig. [TJd) give rough estimates of — T 3 X3 at 
or near T = T sg . We find ~T 3 X3/Xsg varies, as shown 
in Fig. [7p from —T^xz/Xsg ~ 10~ 2 for cubic shapes to 
—T 3 X3/Xsg ~ 10 2 for long thin needles. 

Our results for DID systems with an x = 0.35 site 
occupancy rate can be generalized to smaller values of 
x. As discussed in Ref. Q, any physical quantity / 
satisfies f{x,T) = f(T/x) for x quite smaller than x c , 
the critical concentration above which there is magnetic 



order at low temperature (e.g., x c ~ 0.65 for sc lattices^ 
and 2 ^ x c ~ 0.25 for LiHozYi.^). 
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